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Abstract 

We study the nonlinear hydrodynamics of a 2+1 dimensional charged conformal fluid 
subject to slowly varying external electric and magnetic fields. Following recent work on 
deriving nonlinear hydrodynamics from gravity, we demonstrate how long wavelength per- 
turbations of the AdS dyonic black brane solution of 4D supergravity are governed by 
equations equivalent to fluid dynamics equations in the boundary theory. We investigate 
the implications of S'-duality for our system, and derive restrictions imposed on the trans- 
port coefficients of a generic fluid invariant under the S operation. We also expand on our 
earlier work and determine a new set of previously undetermined transport coefficients for 
the conformal fluid with an AdS gravity dual. Quite surprisingly, we discover that half 
of the transport coefficients allowed by symmetry vanish in the holographic fluid at linear 
order in the hydro dynamic expansion. 
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1. Introduction 



An interesting recent development has been the application of AdS / CFT to the non- 
linear hydrodynamics of strongly coupled field theories [1]] (some early applications include 
0,[J|,|].|J]|||7|,[| ; see also the review [BJ for further references, as well as the review |TTJ for 



earlier work focussed on the linearized regime.) One system of interest, relevant to var- 
ious phenomena including superconductivity, graphene, and the quantum Hall effect, is 
that of a charged 2+1 dimensional conformal fluid evolving under the influence of external 



electromagnetic fields. In a previous paper |Tl[] we applied the method of QTJ] to calculate 



a subset of linear and nonlinear transport coefficients in fluids with an Einstein-Maxwell 
gravity dual. In this paper we extend these results to a complete calculation of all the 
transport coefficients in this sector up to second order in the hydrodynamic expansion. The 
main generalization involves allowing the external electromagnetic field to vary in space 
and time.! Along the way, we study the role of S'-duality in these systems and discover a 
surprising cancelation in the linear order transport coefficients. 

It was noted in that there is a natural SL(2, Z) action on the space of 2+1 



dimensional conformal field theories with a U(l) global symmetry. Generically this action 
is a duality: it maps solutions of one CFT to solutions of some other inequivalent CFT. 
However, in some cases it makes sense to consider CFTs invariant under the S element of 
SL(2,Z). Such CFTs arise naturally in the AdS/CFT context: if a CFT has a gravity 
dual that can be consistently truncated to Einstein-Maxwell gravity, then the boundary S 
operation is mapped to electric-magnetic duality in the bulk, which is a symmetry of the 
Einstein-Maxwell equations of motion. The implications for AdS/CFT linear transport 



properties were discussed extensively in [13,14]. In section 2, we will be interested in 



studying the significance of this invariance for generic S'-invariant CFT hydrodynamics. 
We show how S invariance tightly constrains the hydrodynamics of a generic CFT, so 
that in general, the complete fluid equations of motion at leading order in derivatives are 
completely specified by the equation of state together with one non-negative real function 
and one real number. This is a general statement about S-invariant conformal fluids, 
independent of the AdS/CFT correspondence. 

In sections 3 and 4, we use the dyonic black brane solution of 3 + 1 dimensional 
Einstein-Maxwell gravity to construct a gravity dual for hydrodynamic fluctuations. Fol- 
lowing |I| we proceed order by order in a derivative expansion to arrive at a solution in local 
but not global thermodynamic equilibrium. The starting point is a configuration in global 
equilibrium, with free parameters corresponding to energy and charge densities along with 
the value of a constant background magnetic field. We then allow these parameters to 



3 A natural further generalization is to allow the fluid to live on a curved geometry, as in 
This generalization is discussed in section 2 to first order. 
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vary slowly in spacetime, and also allow for a slowly varying electric field. Associated with 
these varying parameters are numerous transport coefficients expressing the flow of current 
and stress-energy. By solving the Einstein-Maxwell equations in the boundary derivative 
expansion, we determine all these transport coefficients up to second order in derivatives. 

As in our previous work we consider the case of magnetohydrodynamics, such that the 
external magnetic field is nonzero in any Lorentz frame; that is, we assume B 2 > E 2 at 
all points in the fluid. Besides the length scale associated with the temperature, It ~ A, 
there is then a second length scale set by the magnetic field, Ib ~ -?=, where B denotes 

the value of the magnetic field in a Lorentz frame with E = 0. The fluid degrees of freedom 
are assumed to be slowly varying over both of these length scales. Our approach is to be 
distinguished from that of [|l5l , |l4|1 (see also |16|,|17|] ) , which is more appropriate for weak B 
fields (and is also restricted to the linear regime). 

One result of our calculations is that we discover an interesting cancelation in the 
linear transport coefficients for the fluid dual to Einstein-Maxwell gravity. These transport 
coefficients are naturally grouped into two sectors according to their behavior under the 
S operation: the energy diffusion and electric conductivity sector, which was studied in 
T5|,|l4|; and a sector describing charge diffusion and the response to variations in the 



magnetic field. In each of these sectors we find that half of the transport coefficients 
vanish when written in the natural S-covariant basis. This is surprising, inasmuch as we 
demonstrate that nonvanishing values are consistent with all the symmetries and with 
positive divergence of the entropy current. The vanishing results therefore seem particular 
to fluids with an Einstein-Maxwell dual description. Of course, it would be especially 
interesting to find a real world fluid with such vanishing — or at least small — transport 
coefficients, as such a fluid would provide a very promising system for making experimental 
contact with AdS/CFT. 



2. Hydrodynamics of a general ^-invariant conformal fluid 

We are interested in studying the hydrodynamics of a 2+1 dimensional CFT with 
a conserved current coupled to slowly varying external electric and magnetic fields. We 
consider the case where the external field is locally magnetic; that is we impose a constraint 
that \E\ < \B\ everywhere so that there will always exist a boosted frame at each point 
where E = 0. Our fluid dynamical variables are taken to be the energy density, T 00 = e, 
and the charge density, J° = p. In particular, we do not take the fluid velocity to be an 
independent degree of freedom; the reason for this, as discussed extensively in [111] , is that 



in the presence of a background magnetic field momentum is not a conserved quantity. At 
large B and in the fluid dynamical derivative expansion, the fluid momentum density is 
fixed by the equations of motion in terms of e and p, and so introducing additional degrees 
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of freedom to represent it would be a needless complication. 

The equations of motion for our fluid are the conservation equations d v T^ v = F^ v ' J v 
and d^J^ = 0. We solve the equations by working in the neighborhood of an arbitrary 
point in the fluid, and choose a Lorentz frame such that E = at that point. Away from 
the chosen point the electric field will be nonzero but small, in the sense that it will scale 
like a derivative in our fluid dynamical derivative expansion; in particular, its magnitude 
will be of the same hydrodynamic order as derivatives of e, p and B. In this frame it is 
natural to write out the equations of motion as 

= -E % r , 

d^T^ = —pE l + Be ij Jj , 
d p = -diJ 1 , 



d B = e ij d,E 
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with Fqi = Ei, Fij = Beij, and in the last line we have written out the Bianchi identity. In 
this equation we have introduced a notation that we will use throughout, where lowercase 
roman indices refer to the two spatial directions in our fluid. In terms of the quantities in 
(|2.1|) , the SL(2, Z) duality of [[12] corresponds to the following S operation : 



p^B 
B^-p 

E i ^ €ij J 

J 1 -> -e ij E , 



(2.2) 



as well as the T operation: 



p^p + B 
B -> B 

(2.3) 

Ei — > Ei 

r -> ./' • ;'•' /•:, . 

The S operation can be viewed as particle-vortex duality, in that it interchanges a unit 
of electric charge with a unit of magnetic flux. The S and T operations are symmetries 
of the equations of motion (|2.1|) provided that the stress tensor is invariant. Repeated 
application of S and T generates the full SL(2, Z), acting as 



p\ f a b\ ( p \ / Ei \ { a b\ ( Ei 



B ) \ c d J V B J ' \ — €ij J 3 J \c d J V — eij J J 



(2.4) 



4 In the appendix we show how the SL(2, Z) of a boundary CFT can be holographically related 
to the SL(2, Z) of electric-magnetic duality in the bulk for CFTs with a gravity dual. 
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with ad — be = 1. 

In general, S and T map solutions of one CFT to solutions of a different CFT. The 
equations of fluid dynamics are specified not only by (|2.1|) but also by the constitutive 
relations which express the stress tensor and charge current in terms of the fluid variables. 
These relations may or may not be invariant under SL(2, Z). However, it is not hard to 
see that in theories for which the stress tensor is invariant under S, the full equations of 



motion of our hydrodynamic theory will be invariant as well. In fact, as noted in Jl!q] , and 
as we will see in section 3, conformal theories with a gravity dual will be S invariant in 
sectors dual to Einstein-Maxwell gravity. It is these CFTs which we will be considering in 
the remainder of this paper. 

Unlike S, the T of SX(2, Z) will not appear as a symmetry of our equations of motion. 
To appreciate the distinction, note that the action of T is the same as the redefinition 
of charge and current that results from adding a Chern-Simons term to the action of our 
CFT. If we think of theories as being specified by their action, then the T operation relates 
distinct theories. By setting the coefficient of the Chern-Simons term to some fixed value 
(we will take it to be zero), there is then no constraint on the constitutive relations that 
they be invariant under T. By contrast, the S operation relates the theory to itself, and 
we will demand invariance under S. From the bulk gravity point of view, T corresponds 
to shifting the #-angle, while S corresponds to electric- magnetic duality ||12|| . 

To specify the equations of motion for a general S invariant fluid we need to express 
the stress tensor in terms of an S invariant derivative expansions. It turns out that the 
combinations 

M t = tfij =F Be i:j 

3 , J 2.5 
N? = E t T Ji , 

will be particularly useful. These quantities transform under S as: 

(2.6) 



These objects are even under parity, and odd under charge conjugation C. Time reversal 
interchanges the + and — representations. Because the currents and electric field both 
vanish in the hydrodynamic limit, is zeroth order in our derivative expansion while 
N^ 1 is order one. 

There are only two zeroth order S'-invariant scalars in this theory, the energy density 
e in the energy-electric field sector and the quantity 

M = \m+M- % =p 2 + B\ (2.7) 
in the charge-magnetic field sector. 
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2.1. Constitutive equations and transport coefficients 

To completely specify our fluid dynamics we just need to write down constitutive 
relations for and in terms of fluid variables e and p, and the background fields 
B and E$ This is carried out order by order in a derivative expansion, where e, p, B, 
and Ei are all allowed to slowly vary. At zeroth order in derivatives we just have the 
equilibrium fluid, and we can assume without loss of generality that Ei = at lowest order 
in derivatives by locally boosting to an appropriate frame. The equilibrium fluid is thus 
labeled by (e, p, B). 

Given our assumption that B is nonzero, we do not actually need to provide a consti- 
tutive relation for the current. This follows since we can solve for J 1 in the second line of 

(O «* 

J i = -^e ij (d fl T^ + pE j ) . (2.8) 

The current is thus completely determined in terms of other quantities. Note that the 
choice to solve for J; is not invariant under S; the S-dual choice would be to solve for Ei, 
which is allowed provided that p is nonzero. The point is that although our fluid equations 
of motion are assumed to be S-invariant, we are choosing to solve them in a non-manifestly 
S-invariant manner. This asymmetry under S-duality is built in from the assumption that 
B and Ei are regarded as fixed external fields, while p and J 1 are regarded as dynamical 
variables, even though these two sets of quantities are exchanged under S-duality. 

The form of the stress tensor is constrained by symmetries. We demand invar iance 
under charge conjugation, under which and are both odd. Also, we require 
invariance under spatial SO (2) rotations and parity; this requires that i type indices match 
up on both sides of any equation. Note that we do not demand invariance under time 
reversal since, by definition, dissipative fluid dynamics has a preferred direction of time. 
Finally, we impose the scale invariance of our theory, which implies tracelessness of the 
stress tensor, and hence fixes T lJ = \eb~ij at zeroth order in derivatives. Because T 00 = e 
at all orders in the derivative expansion, this leaves only T° l and the symmetric traceless 
part of T 1 - 7 to be determined. 

At first order in the derivative expansion we can define the following four vectors 



We take the fluid variables to be e and p rather than their thermodynamic conjugates T and 
p. This is appropriate, since it is only the former that have an unambiguous meaning away from 
thermal equilibrium. For small deviations from equilibrium one can of course go back and forth 
between the two using the equation of state. 
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sharing the same symmetries as T° l : 



x\ = M+N- 



x f = M±d k M: 



(2.9) 



X 



M- 7 d k M+ 



i ij 

The symmetries then fix T° l to be of the form 

T° l = c x x\ + c 2 x 2 l + c 3 xf + c 4 xt , (2.10) 

where c a = c a (e, M) are, by definition, transport coefficients. We can use dimensional 
analysis to express each of the coefficient functions c a as a power of e times an arbitrary 
function of the dimensionless combination e 4 /M 3 . In the formula (|2.10|) T 0z appears to 



be a function of the current J 1 ; however, it is implicit that J 1 at first order is reexpressed 
in terms of first derivatives of (e, p) through Q2.8| ). Hence T° l is really a function of the 
fluid variables (e, p) and the background fields (B,Ei). There are no contributions to T lJ 
at this order in derivatives. 

2.2. Entropy current 

Studying fluctuations of the dyonic black brane solution will lead to specific functions 
c a , along with their analogs at second order. But before turning to that computation it is 
useful to consider the general constraints on these functions. A set of fluid equations is only 
physically permissible if it is possible to define an entropy current whose divergence is 
positive semi- definite, d^S^ > 0. This expresses the condition that locally entropy should 
be produced and not destroyed. By definition we have S° = s, where s = s(e, M) is 
the entropy density of the equilibrium fluid. To establish consistency we then need find 
the spatial components S l compatible with positive divergence; existence of these spatial 
components will be seen to imply constraints on the transport coefficients c a . As with T° l , 
the symmetries constrain S l at first order in derivatives to be of the form 

S 2 = a\x\ + ot<2%i + ctsxf + a^xj , (2-11) 

with a a = a a (e, M) constrained by dimensional analysis. 

Using the equations of motion, we can compute d^S^ in terms of second order quan- 
tities and squares of first order quantities. Positivity of the divergence for arbitrary con- 
figurations requires that all the second order linear contributions vanish, and that the 
remaining terms form a sum of squares with positive coefficients. 
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Computing the terms linear in second order quantities we find 

[^1 linear = («i - d e s Cl - dMS^jM+m- + (a 2 - d e sc 2 + d M ^jM~ j d i N+ 
+ (a 3 + a 4 - d £ sc 3 - d e sc 4 )M^did k M~ k . 



(2.12) 



Setting this to zero fixes three combinations of the four coefficients a a . Using the parameter 
X to label the family of solutions, we write 



a.\ = d e sci + &ms 
«2 = d e sc 2 - d M s 
a 3 = d e sc 3 + x 
a 4 = d e sc 4 - x ■ 



(2.13) 



With this result in hand, we can proceed to compute the terms quadratic in first order 
quantities. After considerable algebra we find 



+ 

+ 
+ 



d s 

-c x d € s - d M d e s - 

d s 

-c 2 d e s + d M d e s + 



AM 



X i X i 



2 2 

x i x i 



x 

c 3 d M d e s + — + d M X 
X 

c 4 d M d e s - — - d M X 



X i X i 



4 4 

x i x i 



+ [-(ci + c 2 )d 2 s] x\x 



1™2 

i 



+ 



-c 3 d 2 s + d 2 M s - d e x + C!d M d e s + 



9m s 
M 



1 3 

x i x i 



+ [-c 4 d 2 s + d 2 M s + d e x + C!d M d c s] x\x\ 
+ [~c 3 d 2 s - d 2 M s - d e x + c 2 d M d e s] x 2 x\ 

d M s 



+ 



-c 4 d e s - d M s + d e x + c 2 d M d e s - 



M 



2 4 

x i x i 



+ [(c 3 + c 4 )d M d e s] x\x 



3^.4 

i 



(2.14) 



Non-negativity of the divergence of the entropy current now reduces to the condition that 
this quadratic form, viewed as a symmetric matrix, have no negative eigenvalues. 

We can examine this condition separately in various subspaces. In the x\ — xf we find 
the condition ^ 

ci-C2 = -^(2a M ^+^) . (2.15) 
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Next we consider the xf — xf and (xj + xf) — (x| + xf) subspaces, which leads to the 
two conditions 



d M (M X ) 



1 



(c 3 - c 4 )Md M d e s 



d e (M X ) = \d M s - ^(c 3 - c 4 )Md 2 s . 
Equating mixed partials gives an integrability condition 



(2.16) 



d 2 sd M - d M d f sd f 



M(c 3 - c 4 ) 



(2.17) 



This equation determines M(c 3 — c 4 ) up to the addition of an arbitrary function of the 
form f(d e s), which by dimensional analysis takes the form f(d e s) = cj (d e s) 2 . Then, ( |2.15| ) 
fixes Mx in terms of c 3 — c 4 up to an additive constant. 
Using these results the divergence becomes 

dpS* = -- [d 2 s( Cl + pa)] t 2 +^ [d M d £ s(ci + c 2 ) - d 2 s(c 3 + c 4 )] ut-- [d M d e s(cs + c 4 )] u 2 

(2.18) 

where we are writing x} + x 2 = t and xf + xf = u. The determinant of the associated 
matrix is 

1 



16 



d e s(c 3 + c 4 ) + d M d e s(ci + c 2 , 



(2.19) 



Since this is negative semi-definite we must demand that it actually vanishes, which tells 
us that 



c 3 + c 4 



d M d f s 



dh 



(ci + c 2 ) 



(2.20) 



The divergence then finally takes the form: 



Cl + C2f (d 2 s) 2 + (d M d £ s) 2 



y J- _| ry & 



n 2 



d M d e s 



+ (Xj + Xj 



(2.21) 



Since ^ = — T 2 C, where C is the specific heat, we see that d 2 s < and hence c^S 1 ^ > 
provided that c\ + c 2 is non-negative. 

Let us summarize the constraints imposed on an acceptable fluid dynamics. c\ + c 2 is 
allowed to be an arbitrary, non-negative, dimensionally correct, function of (e, M). c\ —c 2 
is fixed by ( |2.15|) ; c 3 + c 4 is fixed by (|2.20|) ; and c 3 — c 4 is fixed by ( [2.171) U P to the addition 



of 



for some number c. Altogether then, given an equation of state there is one 



M(a eS ) 2 ' 

free real function and one free real number labeling the space of allowed fluid equations of 
motion. 

In the above, we considered a fluid living on a flat three dimensional spacetime, but 
one might wonder whether there are any additional constraints to be found by putting the 
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fluid on a curved geometry. It is easy to see that no such additional constraints can arise. 
At this order in the derivative expansion, one can simply repeat the previous computation 
with indices suitably contracted. The one place where a new constraint would arise is if a 
term proportional to the Ricci scalar were to appear in the expression for V A jS ,/x . Since the 
Ricci scalar is not positive definite, positivity of the divergence would impose the constraint 
that its coefficient be zero. However, this coefficient is automatically zero. This can be 
seen by considering the special case of a fluid on S 2 x R. In this case one can obviously find 
static fluid configurations with constant e and p, and these configuration will obviously 
obey V ^S^ = 0. But evaluated on such a configuration, all possible terms at second order 
in derivatives that can appear in the expression for V M .S M vanish, except for the Ricci 
scalar. So the fact that V M -S M = for such a configuration implies that the coefficient of 
the Ricci scalar is zero. We conclude that our fluid equations can be generalized to curved 
spacetime with no new constraints being required. 

2. 3. Black brane equation of state 

In the preceding discussion the equilibrium equation of state of the fluid is left un- 
specified. We now consider the equation of state corresponding to the dyonic black brane 
considered in this paper. The entropy density is given by s = iva 2 where the horizon radius 
a is given by the largest real root of 

or a 

If we use dimensional analysis to write 

d + c 2 = -F(a 4 /M) (2.23) 

ct 

and introduce the numerical constant c, then we find 

3(a 4 + M) 



ci - c 2 = 



AaM 



^=m^m F(aVM) <2 - 24) 

_ (3a 4 + M) (3a 4 - M) 2 
C3 C4_ 4a 2 M + a 6 M C ' 

Any choice of function F and constant c leads to an acceptable fluid dynamics. As we'll 
see later, gravity picks out one particular choice, 

F(aVM) = ?(l + £\ , c = 0. (2.25) 
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corresponding to the transport coefficients 



_ 3(a 4 + M) 
Cl ~ 4aM 
c 2 = 

3a 4 + M 

C3 = 



( 2 - 26 ) 



4a 2 M 
c 4 = . 

It is interesting that gravity yields this simple result in which two of the naturally defined 
transport coefficients vanish. It would be interesting to know if there is some general reason 
for this to be the case since, as emphasized above, all the general consistency requirements 
on the fluid equations can be satisfied even when these coefficients are non-vanishing. 



3. The gravity dual description 



We now turn to our calculation of the transport coefficients for our fluid with an AdS 
gravity dual description. As developed in HT8j[T5[| , the dyonic black brane solution to 4- 
dimensional Einstein-Maxwell gravity provides a dual description of a finite temperature 
2+1 dimensional CFT at nonzero charge density and under the influence of an external 
magnetic field. In this section, we first review the duality dictionary in the case of global 
thermal equilibrium, and then proceed to develop a derivative expansion that will allow 
us to study the duality in the hydrodynamic regime. 



3.1. The dyonic black brane 

The 4-dimensional bulk Einstein-Maxwell action is given by 



S 



d £a/— # 



4 4 



J? tpMN 



2L 2 



(3.1) 



We work in units where L = 1. This action has equations of motion 



Wmn = Rmn + 3gMN — 2F M pFpf + \qmnFpqF p Q — 



Y N = V M F MN = . 



(3.2) 



The dyonic black brane solution to (|3.2| ) is given by 



ds z = Idvdr - U(r)dv 2 + r 2 dx l dx l 
F = -^dr A dv + Bdx 1 A dx 2 



(3.3) 
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where U(r) is the function 

U(r) = r 2 + - | , (3.4) 

and p, e and B are unspecified constants which will later be identified with the charge and 
energy density of our dual fluid as well as an external magnetic field applied to the fluid. 
The black brane has a singularity at r = which is shielded by a horizon at r = a, defined 
by the largest real root of U(a) = 0. We are using Eddington-Finkelstein coordinates 
(r, v , x 1 , x 2 ), which are smooth across this future horizonJl These coordinates do not 
cover the past horizon, and we correspondingly do not demand smoothness there. This is 
physically sensible in that we do not expect solutions of dissipative fluid dynamics to be 
well behaved when extended arbitrarily far back in time. 
The Hawking temperature of the brane is 

T= ^.^±f. (3.5) 
An Ana 3 

e and p are restricted to values such that T > 0. The chemical potential can be read off 
from the Euclidean black hole solution in terms of the the asymptotic value of A v . Recall 
that to obtain a smooth gauge field on the Euclidean section we must choose a gauge such 
that A v vanishes at the horizon (no such requirement exists for the Lorentzian solution 
due to the different topology). This fixes the asymptotic value of A v and gives p = £ . 

Some further conventions: Latin indices M, N, ... run over all four spacetime coordi- 
nates, while Greek indices p, u, ... run over the three coordinates (v, x 1 , x 2 ). Since v plays 
the role of time on the boundary, we will sometimes use v = x . The boundary theory will 
always see a Minkowski metric, j^dx^'dx 1 ' = —(dx ) 2 + dx l dx\ Indices on the boundary 
stress tensor and currents are raised and lowered with this metric. 



3.2. Stress tensor and current 



The action ( |37T| ) should be supplemented with the boundary terms [|19l , |20[1 



bndy 



If*. 2 

-— 2 / a xv— 7 9 2 



4 JdM 



K 1 



d x\/— 7 



(3.6) 



4 JdM 



Here 7 is the boundary metric and 6 = j^O^, where 6^ v = —\{y iJL n u + V^n^) is the 
extrinsic curvature of the boundary, defined in terms of the outward pointing unit normal 
vector n. 



6 These are related to Schwarzschild type coordinates via v = t + r* (r) with = — . 
gives ds 2 = -U(r)dt 2 + + r'dx'dx*. 



This 
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The conformal boundary metric is defined as 7 Mi , = lim^oo -%7 M „. Also, the boundary 
gauge field is defined as linv^oo A M , in a gauge where n M Am = 0. The boundary stress 
tensor and current are then defined as 



5S 



1 



4 JdM 



(3.7) 



Explicitly ||0| 



J* 



T^ V = — [9^ - 9Y V - 27 Ml/ ] 
2 



(3.f 



Implicit in ( p.8[ ) is the large r limit, as well as a projection of T^ v parallel to the boundary 
(since the orthogonal component does not appear in ( |3.7|) .) 

Electromagnetic gauge invariance implies current conservation, 



= . 



(3.9) 



Invariance under diffeomorphisms generated by vector fields tangent to the boundary yields 
the (non) conservation equation 



(3.10) 



Tracelessness of the stress tensor follows from invariance under diffeomorphisms shifting 
the radial location of the boundary 



flV 



. 



(3.11) 



In particular, the latter invariance follows from the absence of logarithmic divergences 
in the bulk action, the presence of which would necessitate adding a non-diff invariant 
counterterm |19 . 

Applied to the solution ( |3~3| ) we find 



(3.12) 



T"" = diag(e, l -e, \e) 
J» = (p, 0,0), 

which demonstrates that e and p are in fact energy and charge densities. 



3.3. Gravitational derivative expansion 

The dyonic black brane supergravity solution of (|3.3|) represents a fluid in global 
thermodynamic equilibrium. We would like to modify this solution in order to account for 
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long wavelength hydrodynamic fluctuations^ so we begin by considering an approximate 
solution which looks locally like the dyonic brane. We start with ( |3.3| ) but allow e, B and 
p to be slowly varying functions of the spacetime coordinates x M : 



d8<® 
F (0) 

where U(r, x^) is given by 

U(r, x^ 1 ) 



2dvdr 
p(x^) 



- U(r, x^)dv 2 + r 2 dx l dx l 
—dr A dv + B(x fl )dx 1 A dx 2 



(3.13) 



2 p 2 (x f ") + B 2 ( X L 
r H ~2 



2e(x' 



(3.14) 



This approximate solution represents a good starting point for two reasons. First, in small 
neighborhoods, it approximates the true dyonic black brane. In the second place, this 
approximate solution approaches an exact solution in the limit of vanishing derivatives 
along the spacetime coordinates. If derivatives are small, then we ought to approach an 
exact solution to the supergravity equations of motion by solving these equations order by 
order in a derivative expansion, and this is exactly how we proceed. 

But first, we deal with some technicalities in defining the gauge potential correspond- 
ing to the field strength in the second line of (|3.13 ). We might consider an that can 
be written in the form 

p(x») 



A (o) 



-dt + \B{x^)e l3 x l A dx 3 + A^(x")dx c 



(3.15) 



where is purely electric, and Ei = d v Af — diA 1 ^ . For constant p, £?, e and A^ this 



form of A^ reproduces the field strength of ( 3.13 ). In addition, this form is generic about 



(o) 



any point since shifts in the origin correspond to constant shifts in A 

But there is also a problem with the expression ( |3.13|) , because generically derivatives 
of B{x^) do not correspond to derivatives of the mag netic field F$ . If we want to match up 
derivatives then we need to correct A^ in ( |3.13| ) order by order in a derivative expansion. 
For example in order to match dF and d 2 F we need to add terms of the form 



-d y By 2 - —d 2 By 3 
4 y 12 y 



dx — 



-d x Bx 2 - —d^Bx 3 
4 12 x 



dy . 



(3.16) 



We note that it is possible to add new terms to this expression order by order in the 
derivative expansion to match up our definition of B in A^ with the magnetic field and 
its derivatives in the boundary. 



Although we use the word "fluctuations", the disturbances are allowed to have large am- 
plitude so long as their wavelength is large. This is the sense in which we are doing nonlinear 
hydro dynamics . 
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We add corrections to (|3.13| ) order by order in a derivative expansion in order to find a 
solution to the equations ( |3.2| ) . As explained in , the equations of motion can be solved 
"tubewise" by working in small neighborhoods near a given location, say x M = 0. 
Each of these tubes then corresponds to a small neighborhood of local thermodynamic 
equilibrium in the boundary fluid. We expand the metric and gauge fields as 

g = g(°\e, p, B, E) + eg^{e, p, B, E) + s 2 g^(e, p, B, E) + 0(e 3 ) 

A = AW (e, p, B, E) + eA™ (e, p, B, E) + e?A™ (e, p, B, E) + 0(e 3 ) , ' ' ' ' 

where g^ and A^ represent the lowest order solution corresponding to the dyonic black 
brane with variable p, e, B, and E. The derivatives of these parameters (and E itself) 
act as sources for the higher order derivative corrections. We have also introduced in 
this equation a parameter e which formally labels the order of a term in our derivative 
expansion. 

( |3.2|) does not admit solutions starting with arbitrary values of B, p, e and E. Instead, 
we find that these parameters must satisfy constraint equations, which we then interpret as 
the fluid equations of motion in the boundary theory. In general, we expect these equations 
to be modified order by order in a derivative expansion. To allow for this, we express the 
energy and charge densities by an expansion, 

e = e (°)(^) +£e ( 1 )( £x ^) + ... i p = p (°)(ex>*) + ep( 1) (ex f *) + --- . (3.18) 

B and E however are specified by Dirichlet boundary conditions, so we wouldn't expect 
for B or E to be corrected order by order in our derivative expansion. Nevertheless, we 
will formally write B and E in a similar fashion, with an eye towards S-duality, 

B = B {0) {ex^)+eB^\ex^) + -- - , E = eE {1) (ex^) + e 2 E {2) (ex») + ■ ■ ■ . (3.19) 

In the standard formulation where E and B are specified by Dirichlet boundary conditions 
E n+1 = B n = for n > 0. Note that this expansion in powers of e is slightly different 
than the labeling used in section 2 where we just counted the number of derivatives acting 
on e and p. If we are working in a small neighborhood around x^ = 0, then it is convenient 
to set e n>0 (0) = p n>0 (0) = E n (0) = 0. 

The zeroth order solution preserves SO(2) rotational symmetry, and this can be used 
to classify the corrections to the metric and gauge fields. We choose a gauge for the metric 
and gauge field 

A r = 0, g rr = 0, g^glZ >0) =0, (3.20) 
and decompose the fluctuations according to their 50(2) representations: 

A {n) = A (n) dv + A {n) dx i ^ ( 3 _ 21 ) 
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with A^ an 50(2) scalar and A^ an 50(2) vector. For the metric we write 

t,(n) 

{ds 2 f n) = —dv 2 - 2h^dvdr + r 2 h^dx i dx i + 2jf>dvdx l + r 2 a ( i fdx i dx j , (3.22) 

In this expansion, and are 50(2) scalars; j\ n ^ is an 50(2) vector; and is an 
50(2) symmetric traceless tensor. 

We impose large r boundary conditions on the n > components given by: 

~ — A (n) ~ - 

r r (3 23) 

jfe(»)~r° a (n) ~- 

These conditions follow from a combination of the asymptotic AdS boundary conditions 
along with the freedom to redefine coordinates as well as the zeroth order solution, as in . 
In addition to these large r boundary conditions, we must also demand that our solution 
be smooth across the future horizon at r = a. For linear perturbations, this condition is 
equivalent to demanding the presence of purely ingoing modes at the future horizon |H 



If we now plug the expansions in ( |3.17| ), ( |3.21| ), ( |3.22| ), into ( |3.2| ) we arrive at the 
dynamical equations for the metric and gauge field corrections. The nth order metric 

(n) 

coefficients are determined by the components of the Einstein equations W^ IN = with 
M, N v . These equations can be organized as 

= -ld r (r 4 cU^) - S$ = 
r 2 (UW rr )^ - Wi n) = d r (™fc (n) ) + d r (d r {r 2 U^)h^) - ^B 2 h™ + 4p^d r A^ - S$ = 

= \rd r (JsdrirjjT^ + - Be^Af - = 

~ \** W ® = d r (~r 2 U^d r a^ - S\f = 

(3.24) 

The source terms denoted by S are constructed from the solution at order n — 1, and so 
are assumed to be known. Similarly, two components of the Maxwell equations yield 

y(n)v = }_Q r ( r 2 drA {n) _ 2p (0) h (n)\ _ y{n) = q 

\ ( 1 \ ( 3 - 25 ) 

= ±d r (u^d r A^ + + B^A ~ Vt ] = 

with source terms V. These equations are sufficient to determine the metric and gauge field 
corrections for arbitrary p, e, B, and E, and are referred to as the dynamical equations. 
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The remaining Y r and W v m are interpreted as constraint equations and will be shown to 
be equivalent to the nonlinear magnetohydrodynamic equations of motion up to second 
order in derivatives. 

Once the metric and gauge field corrections are determined using ( |3.24j ) and ( |3.25|) , 
the current and stress tensor can be found using (a large r limit is implicit): 

r = P 

ji = _ r 2 ^ drA M _ E . 



T 



vv 



(3.26) 

A 



3 ^ („) 



9 A £-~< 



n 



4. Solving the equations 

4-1. Zeroth order solution 

At zeroth order we use the dyonic black brane solution with e = e^ **, p = p(°\ 
B = B(°\ and E = 0. By construction, this is a solution to the Einstein-Maxwell equations. 
The current and stress tensor are 

J" = (p,0,0) 

1 1 (4-1) 
= diag(e, -e, -e) . 



2 2 



4-2. First order solution 
At first order we write 

e ( x ^) = e (°) +£^<V (0) (°) 
p(xn = P {0) + ex^°\0) 
B (0) (^) = S (0) +ex»d^B i0 \0) 
Al (°V) = ^M« (0) (0) . 

The first order sources are built out of d^°\ d^p^°\d^,B^ and and read 

c(i) _ c(i) _ c(i) _ c(i) _ T/(l) _ n 



V 0) + ey«iB(°) _^M- (4-3) 



* r 4 r 4 
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All fluctuations with no sources are set to zero by the boundary conditions. The non-zero 
uations which must be dete 
Integrating Y^ 1 = gives 



fluctuations which must be determined are d r A^ and j± . 



L y " ! (J r A) ' ' + ' C "' ' 

(1) • 



uWdrA™ + + = , (4.4) 



where q is a constant of integration, independent of r. 



Proceeding as in |IlJ] it is easy to obtain 



(J) _ Mr.jj^+rdjM-j+r 2 ^ 



(4.5) 



d r A\ L \r) = 



! C/(°)(r) 



We fix the integration constants by imposing (|3.23|) and demanding that the fluctuations 



be smooth across the future horizon at r = a. This fixes ac- = oo and 

=-±(l+ ) c« —*-(l+ ) 9 M- (4 6) 

Pl a(°) V a(°)U"(°)(a(°)) J 1 a (o) 2 V + a(°)£/"(°)(a( )) J 3 lJ ' 1 J 

With this choice of integration constants we find that has the large r behavior 

> (i) (r) = ^L +0( ^ ) . (4.7) 
To find the current J 1 we need 



,2 F ir (i) = _Af 7j f) _ jgW _ ^ _ C/Wa^W . (4i 



At large r, we use this expression and ( |3.26|) to find 

ji = c (d _ S (D ) (4 . 9) 

which tells us that 

Cl (1) = N~ . (4.10) 
The stress tensor can be found at this order using ( |3.26| ) 

4 13 3 (4.11) 
[T ij ] st = . 
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where the symmetric traceless part of a matrix is defined according to 



[M, 



1st 



, yJ — - (Mij + Mji — SijMkk) 



(4.12) 



If we now combine (|4.6| ) and ( fl.ll|) with our definition of the transport coefficients 
( p.lO| ), T° l = c\x\ + C2xf + C3xf + c^xf, then it is easy to derive (|2.26| ) 



ci 



3(a 4 + M) 
AaM 
c 2 = 

3a 4 + M 



(4.13) 



C3 



4a 2 M 



c 4 = 



As advertised, we find the surprising result that two of the transport coefficients vanish. 

We interpret the remaining W v m an d Y r equations as constraints on the allowed 
values of the fluid variables. If we use our solutions for and d r A^ found in equations 
) and Ql5|), then it can be shown that these equations reduce to 



(4.14) 



at order e in our derivative expansion. In particular these constraint equations take the 



form 



d vP W = &,e<°> = 



J 1 



€jjdn€ — —6ioE 



.(1) 



1 

1 

-e 



P_ 



B(°) 



ij 



(4.15) 



4-3. Second order solution 

e, £?, p and A^(x^) are now given by expanding to order £ 2 , 



e(x h 



p(x^ 



= e (0) (0) +£x^a M e (0) (0) + -e 2 x^x v d^d v t {0 \Q) + £ V^eW(0) 
= p (0) (0) + £x^p (0) (0) + -e 2 x fl x iy d^d J/ p w (0) + e 2 x fl d^p {1) (0) 
S(x M ) = B(°)(0) + 6x^8^3^(0) + -e 2 x^x 1/ d^d u B^\o)+e 2 x ll d ll B^\o) 
Al (x») = ex^d^Al (°)(0) + -e 2 x^x v d^A^ (°\o) + e 2 x»d fI A% «(0) . 



(4.16) 
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The second order sources work out to be 



;(!)\2 



5{ 2) = o 



o(2) 



(1) P .(1W 



(2) 



. 



(4.17) 



It is not hard to work out the solutions to these equations following |Tl] 



d r A&(r 

fe (2)( r 
■(2)/ 



(2)/ 



-2 fdr" (r"d r Af\r")\ 

J oo T J oo 

i r 

—~ / dr'r' 2 Xi(r') 
\r [ dr'X 2 (r') 

^ Joo 



-U (0) (r)M± I dr 



*(2) , 4 (2) 



C7(°)(r 



A2 



(4.18) 



where 



X 1 (r) = -p rh W 

X 2 (r) = d r (d r (r 2 U^)h^) - ^S 2 /i (2) + Ap^d r A^ - 5, 



? (a) 
'(fc) 



(4.19) 



The scalar sector does not contribute to the currents and we will not study the scalar 
sector solutions in any greater detail. 

In the vector sector, d r A\ is determined by the analog of ( |1.4| ), 



2 iJ J J 



-c„- 



(2) 



(4.20) 
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Imposing regularity at the horizon tells us that 



(2) 



a 



(0) 



1 + 



U'^(a^) \ (2) 



a( )C/"( )(a(°); 



(4.21) 



(2") 

01 ' and cf' now contribute to the stress tensor and current precisely as in (|4.9|) and 



.(2) 



Also, the stress tensor conservation equation now fixes i 



.(2) 



1 



(2) 



N 



(2)- 



(4.22) 



The vector components of the stress tensor and current at this order take the same form 
as at first order. 

In the tensor sector, the solution in ( |4.18| ) is straightforward, and leads to 



1 r°° 

[T ij ] st =» drS%\r). 

2 Ja(0) 



(4.23) 



In order to write down the transport coefficients at this order we want to make use 
of the symmetries to organize the stress tensor according to the S'-duality representations. 
We use the operators K^j = M^M^- with eigenvalue —1 under S. Below we list all of the 
terms which appear in the symmetric traceless part of the stress tensor at second order 
in derivatives, there are of course other terms consistent with the symmetries, but we list 
only the ones which appear in the fluid from our gravity dual. 



[Tij yt =b NN [ K + N - N -]« + b NN [N+N-f 



+ b NM 



K+NZd,M 



3 1 



si 



NM 



+ 6_+ 



n st 



Nr dl M+ 



NM 



-i st 



n st 



d m M±,dlM. 



n st 



+ b N [M+W-] st + b M M+d^M- 



st 



3 1 



(4.24) 



To read off the expressions for the transport coefficients we use ( f4.23|) and the expressions 
above. We begin with the two linear transport coefficients 

2M (4.25) 
b M = " 



2M 



In many instances in this paper, we have used simply TV* and to refer to N- 1S>± and 



M^ ^; the difference is higher order in the derivative expansion and so we have not been overly 



careful to distinguish between the two. Strictly speaking, we have = My +• ■ -, and 

it is implicit that we are to take the lowest order when these quantities are present in derivative 
expansions. 
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As written in ( 4.24 ) it may appear that these are only half the coefficients consistent with 
the symmetries, but there are only two linear transport coefficients present at this order 
in an arbitrary S'-invariant fluid. We emphasize that there is no mysterious cancelation in 
([125]), unlike what we found at first order. 

The nonlinear coefficients are more difficult to obtain. The calculation is straightfor- 
ward but not especially illuminating, so we simply state our result: 

1 r°° ( 9 9 9 

■NN L , 1 / 7 / z • •/ , z - • , 1 •/ „•/ z 



a 



b — = 2nMT + 2 / dr [~r jNjN + ^ 3n3n + * 3n3n ~ ^ nAn 



b MN = 3(a 4 + M) f°° ^ f j M j' N _ JjvJm , Jn3m i ., ., _ JnAm + 3m A N 

AtKiM 2 T J a \ r r r 2 2 N M r 2 

— = 8rrM 2 T + 2 J a \~r 3M3M ^ 3M3M + * 3m3m ~ ^2 JmAm ) 



6 - = + J. * A » 



X 



■MM 1 , / 



*" 87ra 2 T 7 a r 

(4.26) 

where the prime notation indicates a derivative with respect to r and we have written the 
vector correction to the metric and gauge field at first order (|4.5| ) as 



with 



j( 1 \r)=j N (r)M+N- + j M (r)M+d k M7 k 
d r Af\r) = A N (r)N~ + A^djM^ , 



U'{a) 



(4.27) 



3n{t) = / dr' 



a Joo U*{r) 

nn/ \ r r 1 _1_ 2U '( a ) a 

]M(r) = !£W \ d /. + W ~ > 



a- 



X 



U 2 (r) (4.28) 



A / N P 2 +B 2 . 1 



t/(r)r 2 JJV C7(r) 

4 r ^ P 2 + g2 - 1 



U(r)r 2 U(r)r 

Unlike the linear second order transport coefficients, there appear to be fewer coefficients 
present in the second order nonlinear sector than are predicted by symmetry. Coefficients 
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of terms like [K+N+N+} st and [K~ k N+ N+] st are seen to be zero, and it is not hard to 
trace through our calculation to see that the cancelations at first order are responsible for 
the vanishing of these terms. While it is possible that some (or even all) of these missing 
coefficients might be explained through a rigorous analysis of the entropy current at second 
order, we do not yet have an understanding of why some representations of the ^-duality 
appear to be favored at this order. 

In terms of these currents, it can be shown that the remaining constraint equations 
at this order reduce to the equations of motion in our fluid (|2.1|) . 

5. Conclusion 

Let us review what has been achieved. We first considered the fluid dynamics of a 
general S-invariant fluid at first order in the derivative expansion, independent of gravity or 
the AdS/CFT correspondence. Even at this first nontrivial order there are a large number 
of transport coefficients, since we allow for arbitrarily varying energy and charge densities, 
and arbitrarily varying background electromagnetic fields. By using the constraints of 
symmetry and positive entropy production, we found the most general form of the transport 
coefficients, and found that they could be expressed in terms of an arbitrary real function, 
whose argument is the single ^-invariant dimensionless combination of the fluid variables, 
along with one real constant. This result implies many nontrivial relations among the 
various transport coefficients. In principle, it would be possible to extend this analysis 
out to second or higher order in the derivative expansion, although the number of terms 
proliferates rapidly. 

We then turned to the gravitational description of fluid dynamics in terms of black 
branes in an asymptotically AdS4 geometry. By solving the Einstein-Maxwell equations 
order by order in a boundary derivative expansion we were able to compute all transport 
coefficients up to second order. At first order our results were in agreement with those 
expected from our general fluid analysis, and we obtained the specific forms of the free 
function and constant that appeared in that analysis. Interestingly, this yielded vanishing 
values for two of the four transport coefficients, working in a basis natural under S-duality. 
At second order we obtained new results for both linear and nonlinear transport coefficients. 

S-invariance is a property shared by any fluid that has a holographic description 
in terms of four dimensional Einstein-Maxwell theory. One can hardly resist wondering 
whether such a fluid exists in nature. A signal that one had found such a fluid would be 
to verify the relations among transport coefficients that we derived in section 2. An even 
more remarkable result would be to find a fluid with vanishing (or small) values for the 
two transport coefficients which vanished in our gravitational computation. This would 
be a smoking gun for holography, since this vanishing does not appear to be fixed by any 
obvious symmetry or consistency condition. 
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Appendix A. S and T operations: bulk versus boundary 



In section 2, we referred to the observation of |H| that 4-dimensional electric- magnetic 



SL(2,Z) duality in the bulk is mapped to the operations we defined in (|2.2| ) and 
The purpose of this appendix is to establish this. 

We begin with the S-operation and consider the most general field strength allowed 
by our boundary conditions and choice of gauge. Using ( |3.13|) , (|3.15| ), and fl3.17j) , we are 
restricted to considering an F of the form 

F = A^dr A dv + Bdx 1 A dx 2 + Edv Adx 1 + d M A^\r)dx M A dx N , (A.l) 

n 

where all of the elements above are assumed to be functions of the boundary coordinates 
x M . Using ( |3.26| ) and the boundary conditions we can write this as 

£+0(^))drAA; 
+ (b + 0(-)) dx i Adx J 
+ ( Ei + 0(-)j dv Adx 1 
+ (-^*±^. + +0(\)) drAdx* 

Given the allowed form for the metric, as outlined in section 3.3, we compute the dual field 
strength 

*F = ( 4 + 0(4) J drAdv 



+ (-p + 0{-)} dx 1 Adx 2 

. 1 \ . (A-) 

+ {-eijJ J +0{-)\dvAdx l 

( fiJ _|_ J 3 1 \ 

+ ( e v — ~2 — + drAdxl ■ 

Comparison of (|A.2|) and ( |A.3| ) demonstrates that electric-magnetic duality in the bulk 
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exchanges the boundary CFT parameters according to (2.2) 



p^B 
B^-p 



J 1 -e ij Ej 



(A.4) 



Turning now to the T operation, we consider adding to the action ( |3.1| ) a 9 term of 
the form 

Se = J~i [ e M NOpF MN F OP . (A.5) 
According to ( |3.7| ) this term modifies our definition of the current to 

= r 2 ( F r » + ^-e r ^ p F v A , (A.6) 



e 

where a large r limit is implicit. In the bulk, the T operation corresponds to 9 — > 9 + 2tt. 



In the boundary theory, this corresponds to (|2.3| ) 



p^p + B 
B -> B 

(A.7) 

J 1 -> ./''•- <''•'/•;,• . 

This establishes the relation between the £ and T operations in the bulk and boundary. 
It also shows the different status of the two operations: S corresponds to a symmetry of 
the bulk equations of motion, while T corresponds to a change in the action of the theory. 
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